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We analyze the effects of hyperfine interactions on coherent control experiments in triple quantum 
dots. By exploiting Hamiltonian symmetries and the SU(3) structure of the triple-dot system under 
pseudo-exchange and longitudinal hyperfine couplings, we provide analytic formulae for the hyperfine 
decay of triple-dot Rabi and dephasing experiments. 
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I. INTRODUCTION 

Recent experiments in triple quantum dot flEI show in- 
creased coherent control abilities over previous experi- 
ments in double dots^^. In both cases, one factor lim- 
iting coherence is the hyperfine interaction to substrate 
nuclear spins. In this brief report, we derive simple ex- 
pressions for the form of hyperfine-induced decay in sim- 
ple triple-dot experiments, and indicate how these ex- 
periments may be used to elucidate the total volume of 
the quantum dot wavef unctions. The formalism we in- 
troduce may also impact the development of composite 
control sequences for mitigating hyperfine-induced pulse 
errors. 

In our analysis, we treat the nuclei as providing clas- 
sical random fields. In an applied magnetic field suffi- 
ciently large to energetically forbid electron-nuclear flip- 
flops, the hyperfine Hamiltonian may be well approxi- 
mated as Hhf — J2j BjSj, where Sj is the z-component 
of the spin operator for spin j, with S = 1/2. Here 
z refers to the direction of the applied magnetic field. 
The hyperfine fields Bj are summations over all nuclear 
spins within each quantum dot volume, with appropriate 
coupling constants^'. These sums are treated as clas- 
sical random variables drawn from a zero-mean Gaus- 
sian distribution of variance (Bj) — en. The average 
value of <7j across three dots, which we notate cr HF , pro- 
vides our scaling throughout this analysis. In r eal time 
units, o-jjp 1 ~ T 2 ~ 10 ns f° r GaAs-based dotiP^ and 
about 300 ns for Si-based dots 8 . Nuclear fluctuations are 
driven by nuclear dipole-dipole kinetics which are very 
slow in comparison to typical control and measurement 
timescaleiP. 

Coherent control of the spins is achieved by pseudo- 
exchange between dots 1 and 2, or between dots 2 and 3. 
Exchange operators are written in terms of spin operators 
as Ejk = Sj ■ Sfc. The exchange terms of the Hamiltonian 
may be modulated via voltage pulses. 

Naively, the analysis of hyperfine-averaged dynamics 
under various exchange pulses may be performed by con- 
structing the full Hamiltonian in an 8-dimensional basis 
and then averaging over the 3 independent random vari- 
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FIG. 1. (a) Schematic of energy levels of triple-dot system, 
above a sketch of triple-dot energy potentials. The black dot 
and black dashed arrow indicate the initialization and puls- 
ing of a triple-dot Rabi or dephasing experiment, (b) SU(3) 
structure of triple dot system. 



II. SU(3) DESCRIPTION OF TRIPLE-DOT 

The first important symmetry of both -ffjjF an d the 
exchange operators is that they all conserve the z- 
projection of the total angular momentum of the three 
spins. It follows that the Hamiltonian may be block- 
diagonalized into four blocks corresponding to quantum 
number m z = ±3/2, ±1/2. Since current experiments 
initialize into a spin-singlet state, the two m z = ±3/2 
states are never occupied, and hence we will discuss them 
no further. The remaining six states are divided into two 
manifolds, which we notate ft and ft for m z = 1/2 and 
m z = —1/2, respectively. We refer to this degree of free- 
dom as a gauge qubit. 

A convenient set of E12 eigenstates in the ft manifold 
are 
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The states are notated according to logical states of a 
qubit in a decoherence free subsystem (DFS) 10 * 11 ! The 



ables Bj. Employing three symmetries, we reduce this labels -D±i/2)> etc. are those of Refs. [T]andHl provided 



problem to a one-dimensional integral with a good an- 
alytic approximation; the resulting expressions are then 
particularly useful for curve-fitting, allowing the estima- 
tion of quantum dot volumes via hyperfine parameters. 



here for comparison. The ft manifold may be found by 
flipping all three spins in the ft manifold. This choice of 
states is physically motivated by the ability to initialize 
into the ground state of the exchange operator E12 by 
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preparing a ground state in the (2,0,1) charge configura- 
tion, as shown in Fig. [lja). This ground state is chosen 
as logical 0. The state labeled Q has the critical feature 
of being fully symmetric and therefore equally unaffected 
by exchange on any pair of dots; this is therefore not use- 
ful for coherent control and is referred to as a "leaked" 
state. The other logical qubit state is the only remaining 
orthogonal state within the given gauge. 

Figure [lja) sketches energy levels in only one of the 
gauge manifolds, for simplicity; the other one possesses 
the same structure. We therefore drop the ft,])- distinc- 
tion from the state labels. More complete diagrams may 
be found in Refs. [Hand El On the left of this figure, cor- 
responding to negative detuning parameter e, the three 
dots have the (2, 0, 1) charge configuration, and as e is 
increased, hybridization with the (1, 1, 1) charge configu- 
ration results in pseudo-exchange of magnitude J z , in- 
dicated as the energy difference between the red and 
blue lines, corresponding respectively to the |0) state |1) 
states. The states |0) , |1) , and \Q) are energy eigenstates 
only far to the left of this energy diagram. The subscript 
z is chosen for J z since the Hamiltonian term J Z E\2 splits 
the qubit states, and is therefore logically analogous to 
the Pauli operator <j z on the logical qubit. Far to the 
right of the diagram, pseudo-exchange between dots 2 
and 3 has magnitude J n . 

The choice of n subscript on J n is made evident by 
considering the group structure for a single gauge. All 
operations of exchange control and hyperfine drift within 
one gauge are described by the group SU(3). Using the 
states in Eq. (JlJ in the order given, a natural selection 
of generators are the Gell-Mann matrices in each basis, 
traditionally labeled Xj for j = 1 ... 8. These matrices 
may be interpreted as raising and lowering through a 
triangle of basis states as diagrammed in Fig. [ljb) . The 
exchange operators are independent of gauge, and may 
be written in terms of the SU(3) generators as 



written 



E\2 — 

E23 = 



2 v / 3 



2\/3 



A3 cos <j) + Ai sin < 
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where <fr — 120°. This representation makes evident the 
effect of each term for the DFS qubit Bloch sphere, de- 
fined by the projections to Ai,A2,A3. The As in each 
term commutes with these generators; it provides a rel- 
ative phase to the |Q) state. The E12 term provides 
rotations about the A3 or i axis, and hence we ascribe 
it the magnitude J z . The £"23 term provides rotations 
about the h = [sin(/>,O,cos0] axis, and hence we ascribe 
it the magnitude J„. The exchange Hamiltonian H ex = 
JzE\2 + J n E23 is diagonalized by the unitary ^2(77) = 
exp(— 177X2/2), where 77 = tan -1 [J n sin <j>/{ J n cos 4> + J z )], 
providing the eigenstates |A ±1 / 2 ) and |A' ±1 ^ 2 ) of Ref.[T] 

In our SU(3) picture, the hyperfine Hamiltonian is 



H m = ±U 7 tf) 



A, 



A 1 + -^A 8 )l4(/3), (4) 



where U 7 (f3) = exp(-iX 7 /3/2) for = n - tan" 1 V8 « 
109.5°. The hyperfine differences are A 12 = B\ —B 2 , and 
Aj9 = B3 — (B\+ B2) /2. We will also use the linear com- 
binations A23 and A53, defined via permutation of the 
subscripted dot labels. These hyperfine field differences 
have variances and covariance 
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The action of these hyperfine terms can be assessed using 
Fig.fijb). First, U 7 {/3) mixes the |1 ft) and \Q ft) states 
into jT^ t) and (T+ 1 |), where \T$) is the J z = m pro- 
jection of the J = 1 triplet for spins 1 and 2. These 
states are degenerate with respect to i?i 2 , and hence 
[A7, -E12] = 0, which provides our second key symmetry 
allowing analytic diagonalization. Following this fixed 
SU(3) rotation, the hyperfine Hamiltonian is highly rem- 
iniscent of the double-dot singlet-triplet case: A 12 causes 
only x-axis rotations of the double-dot singlet-triplet 
qubit, while A^ provides a relative phase to the polar- 
ized triplet state. 

The final symmetry we employ is that the exchange 
interaction is independent of the two gauge manifolds, 
ft and JJ-, while Eq. Q exactly flips sign between 
them. We briefly note that this points to hyperfine 
decoupling schemes in which the total electron spin is 
rapidly flipped, perhaps with electron-spin-resonancc- 
based 7r-pulses; such schemes would have no direct effect 
on exchange interactions. For the purposes of our present 
analysis, if an experiment averages over A12 and Ay^, 
and if the probability distributions for these variables are 
symmetric about zero, then the average dynamics will be 
the same in each manifold. This allows us to drop the 
f|\4J. distinction and work with a single SU(3) system. 

There is one further energy term in the system which 
we have not addressed, that of a global magnetic field, 
either experimentally applied or due to average hyper- 
fine terms. A global field splits the energy of the gauge 
manifolds, but plays no further role; it is of course the 
invariance to this field and its fluctuati ons to which the 
nomenclature "decoherence-free" applies^LUJ 



III. RABI AND DEPHASING EXPERIMENTS 

The algebra we have introduced enables an analytic de- 
scription of two of the simplest experiments which may 
be performed in a pulsed triple-dot system. In both ex- 
periments, the singlet |0) is prepared in dots 1 and 2 [see 
Fig. [lja)] with the gauge state typically random. (High 
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magnetic-field experiments may allow for initialization 
into a single gauge state). For a Rabi-like experiment, 
a bias sweep then quickly reduces the magnitude of J z 
while increasing the magnitude of J n . Free evolution with 
J„ > then occurs for time t, after which J z is quickly 
ramped up as J n is quickly ramped down. The proba- 
bility of returning to the initial singlet state is observed 
via Pauli blockade on the (1,1,1) / (2,0,l)charge transition 
and subsequent charge-state detection^. A dephasing 
experiment in a triple dot is the same experiment, except 
J n is chosen to be a very small value, optimally zero. 

In Ref. [TJ the triple-dot Hamiltonian is diagonal- 
ized for nonzero J„ and J z , but with vanishing hyper- 
fine. Here, we diagonalize for a single nonvanishing ex- 
change term and nonvanishing hyperfine. For the ex- 
periments described, the Hamiltonian to be diagonal- 
ized is H = H HF + J n E23- We begin by diagonaliz- 
ing E23 via the rotation U2(4') = ex P( — *<M2/2), since 
Ul(4i)E23U 2 (0) = Ei2- In general, this rotation may be 
considered to permute the dot labels; the rotated hyper- 
fine Hamiltonian may therefore be taken as Eq. Q with 
the {1, 2} subscripts permuted to {2, 3}. Stated another 
way, the problem of initializing as a singlet in dots 1 and 2 
and pulsing J n is equivalent to the problem of initializing 
as a singlet in dots 2 and 3 and pulsing J z , except that 
we must swap the labels of dots 1 and 3 in the hyperfine 
parameters. As we have seen, we may then apply Ut((3) 
to block-diagonalize H , and we are left with diagonaliza- 
tion of a single-qubit Hamiltonian. A final A2 rotation by 
angle 9 = — tan -1 A 23 /J n completes the diagonalization: 
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We then calculate the probability Po{t) = 
(|(0| e-^ HKF+JnE23)t |0)| 2 ) by averaging A 23 and A23 over 
their joint Gaussian distribution; the integral over A23 
is trivial. The final result may be expressed in terms of 
two one-dimensional integrals over A23 as 



Po(t) = 
2 4 V ' 4 



The first integral is 
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FIG. 2. Numeric example of Po(t) for three values of J n . The 
example values for the hyperfine deviations for the three dots 
used for this plot are Oj = (j/2)a HF . The black curve shows 
numeric integration of Eq. pj, featuring no approximation 
relative to full 8x8 diagonalization and three-dimensional 
averaging. T he blue dashed line shows the high-J n approxi- 
mation [Eqs. | 14]|T6 l]. The red line (short dashes) shows the 
low-Jn approximation [Eq. ( 18 1] . For each J n a Bloch sphere 
is shown, in which the three cardinal directions are the expec- 
tation value of the first three Gell-Mann matrices Ai, A2, and 
A3. Hence the surface is the space of DFS qubit states, while 
the interior corresponds to leakage. The 8 colored trajecto- 
ries represent the unitary evolution under 8 different hyperfine 
fields, sampled from the 3 independent Gaussian distributions 
of the 3 dots. Po(t) is the average projection of these trajec- 
tories onto the z axis, from which all trajectories begin. 
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where u = 023*, x = A23/CT23, and y = J n j<J2z- In terms 
of the same variables, the second integral is 
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v^ 2 + y 2 

x 2 /2+i(y/x 2 +v 2 -y)u/2^ 



where w — C2 3 /cr| 3 . These integrals resemble those oc- 
curring in the analysis of double-dot experiments^. Ex- 
amples are plotted as the black lines in Fig. [2] 

Fortunately, the Rabi and dephasing experiments de- 
scribed above function in one of two limits. The Rabi 
experiment operates in the high y — J n /<7 2 3 limit. For 
y ->• 00, h -> 1 and I 2 {t) -> exp[-(C 2 3V (7 23) 2 /2] result- 
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ing in the simple expression 



where 



Uxn m = 5 + cos(^t) + I + cmM,-^ 
j„->oo 8 8 4 

(12) 

In contrast to double-dot Rabi experiments^, this ex- 
pression features an offset and oscillatory term decaying 
with effective T 2 * of V2/a^, independent of J n . The first 
order corrections to the integrals are expressed in terms 
of the width function 



(13) 



as 



1 - y/A(t,J n ,0) cos{J„i + (1/2) cos- 1 A(t, J„,0)} 



(14) 



where F(a;) = v / 27rxexpx 2 erfc(a;), and 



$(t) = - cos- 1 A(t, 2J n , 0) - cos- 1 A(t, 2J„, C 23 /4j) 
-(C 23 /a 2 z ) 2 [l-A 2 (t,2J n ,0)(t)}- (16) 



These functions describe an additional long-tailed decay 
with timescale Jn/<7 23 and associated phase modul ation, 
similar to Rabi experiments in the double-dot system^M. 
These long tails are often obscured by the effects of charge 
noised 

The dephasing experiment operates in the low J n /<J23 
limit. In the limit J n = 0, 



-(<Ti 2 t) 2 /2 
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-A 2 (t.2J„,0)C| 3t 2 /2 CT l3{ cos [ $ ( i )] + cos[J< + ( 15 ) 



exactly the same as in the double-dot cas d 4 * 5 * 9 * since the 
third dot plays no role. The situation is quite different 
from the double-dot case at finite J„; a low-order approx- 
imation to the integrals Ii{t) and / 2 (i) yields 
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This approximation as well as the high- J approximation 
are compared to exact integrals in Fig. [2] 



IV. DISCUSSION 

Several intuitive observations may be extracted from 
our analytic results. First, if one seeks to measure the 
volume of electronic wavefunctions in a triple-dot sys- 
tem via hyperfine effects, one needs enough information 
to explicitly obtain <7i,<r 2 , and cr 3 , in addition to the 
material-dependent hyperfine coupling constants and iso- 
topic content. Our expressions show that the hypcrfinc- 
decay of Rabi oscillations provides a^, while the dephas- 
ing experiment provides ci 2 . These two decay periods 
are sufficient to extract erf explicitly. If signal-to-noise 
and charge-noise concerns allow it, er 23 may in princi- 
ple be extracted from either experiment via curve-fitting, 



although a more reliable measure is to physically pre- 
pare the initial singlet via the (1,0,2) charge state and 
pulse J z , providing aj^ and cr 23 via the dominant decay 
in Rabi and dephasing experiments. Note that similar 
experiments using singlet-triplet control in double-dots 
do not allow extraction of hyperfine parameter s of i ndi- 
vidual dots, only the variance of their differenc e 1 5 * 8 * 9 ^. 

Another important difference between triple-dot Rabi 
and double-dot Rabi is the ability to distinguish noise in 
the exchange parameter, J, typically caused by charge 
noise, from that of hyperfine noise. Equation 12 shows 
two oscillatory terms of comparable magnitude, one 
showing hyperfine decay and one not. Fluctuations of J 
(perhaps due to detuning fluctuations) in averaged exper- 
iments would result in decay of both these terms, but hy- 
perfine decay is only active on one of them, allowing one 
to distinguish the different effects within a single dataset. 
Yet further information about the noise characteristics is 
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available from more complex pulsed experiments, such 
as triple-dot composite pulsing and spin-echo-like exper- 
iments. The SU(3)-based analysis we provide here may 
assist the analysis of hyperfine effects in these more com- 
plex experiments, as well. 
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